function identityC=computeIdentityC(Hp, C)
[rC cC] = size(C);
identityC = zeros(Hp * rC, Hp * cC);

for i=1:Hp
    [rStartPointer rEndPointer] = getPointers(i, rC);
    [cStartPointer cEndPointer] = getPointers(i, cC);
    
    identityC(rStartPointer:rEndPointer, cStartPointer:cEndPointer) = C;
end
end